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ABSTRACT 

We present Arecibo L-band Feed Array 21-cm observations of a sub-complex of HVCs at the tip of 
the Anti-Center Complex. These observations show morphological details that point to interaction 
with the ambient halo medium and differential drag within the cloud sub-complex. We develop a new 
technique for measuring cloud distances, which relies upon these observed morphological and kinematic 
characteristics, and show that it is consistent with Ha distances. These results are consistent with 
distances to HVCs and halo densities derived from models in which HVCs are formed from cooling 
halo gas. 

Subject headings: ISM: kinematics and dynamics, ISM: clouds, Galaxy: halo, radio lines: ISM 



1. INTRODUCTION 

High-velocity clouds (HVCs) have been detected all 
over the Galactic sky in the 21-cm hyperfine transition 
of neutral hydrogen (for an in-depth review see the book 
edited by Ivan Woerden et aHl200l . HVCs are clouds 
with typical velocities between —300 and 300 km s _1 in 
the Galactic Standard of Rest (GSR) that cannot be ex- 
plained by Galactic rotation. Distances to HVCs, and 
therefore their physical scope, are uncertain and con- 
tentious, but recently constraints have been put on the 
distances to some clouds. The rare direct absorption 
line measurements of halo stars typically provide distance 
lower limit s of 2 kpc and include a few upp er limits of 10 
kpc (e.g. iThom et al]l200fl lWakkerll200ll for a review). 
Successful measurements of Ha flux from HVCs, which 
indicate the ionization of the HVCs by escaping Galactic 
radiation, suggest that many HVCs are within 40 kpc, 
though the inferred distances have large uncertaintie s 
^Bland-Hawthorn fc PutmanlfeOOlUPutman et alJl205^ . 
It is crucial that we understand the true physical size of 
these HVCs; they have been invoked to provide the ac- 
creti ng material that fuels sta r formation in our Galaxy 
(e.g. iR.ocha-Pinto et aLll2000fi and may indeed be a sig- 
nature of the primary way that baryons cool to form 
galaxies like the Milky Way. HVCs may also lend in- 
sight into the Galactic formation pro cess as a tracer of 
the e volution of the Galactic halo (e.g. iMaller fc Bullockl 
2004). 

HVCs are thought to move through the Galactic halo. 
Observations of Ovi absorption on sightlines to HVCs 
have shown that there is collisional ionization in the 
HVCs, p resumably from inter action with the ambient 
medium ijSembach et aljr2003|) . Also, many HVCs are 
seen to have 'head-tail' structures, wherein the cloud 
has a colder high column density core, surrounded by a 
warmer HI shell which trails off in one direction from the 
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cloud with an associated velocity gradient, a s if the cloud 
were be ing ablated b y the medium (e.g. IBriins et alJ 
HXXJ. IWolnre et all lll99l (hereafter, W95) showed 
that the two-phase structure of some HVCs implies a 
significant pressure in the ionized halo that declines 
with height above th e Galactic disk. Subsequently, 
iSternberg et alJ l)200 2) showed that this two-phase struc- 
ture implies that HVCs reside in the Galactic halo and 
not the Local Group halo. 

We have obtained an image in the 21-cm line of a 
small sub-complex of high-velocity clouds (HVCs) lo- 
cated at the tip of the Anti-Center Complex with the 
Arecibo telescope 5 . Arecibo paired with the Arecibo 
L-Band Feed Array (ALFA), a large filled-apcrture tele- 
scope with 7-fold multiplexing, affo rds extremely high 
sensit ivity and dramatic resolution (jStanimirovic et all 
2006). This combination allows us to study the mor- 
phology of large clouds and their interfaces with ionized 
ISM phases in fine detail as never before. The data sug- 
gest interaction between the halo gaseous medium and 
the HVC itself. We claim that by modeling this inter- 
action we can gain insight into the physical attributes 
of the cloud and the medium through which it moves. 
Though dr ag has been studied befor e in the context of 
HVCs (e.g. IBeniarnin DaulvllTflfll . it is our intent to 
model differential drag within this cloud group, which 
is immediately evident in the morphology and velocity 
structure. In particular, these calculations can give us 
another handle on the distance to the cloud in question, 
as well as the density of the ambient medium. 

2. OBSERVATIONS 

We observed a region of sky 16x7 square degrees cen- 
tered on a = 2 h 15 m , S = +9°30 m in May and June of 
2005 with Arecibo's ALFA receiver and the GALSPECT 
spectrometer. ALFA is a new 7-element focal-plane ar- 
ray primarily designed for 21-cm observations. GAL- 
SPECT is a special-purpose spectrometer for Galactic 
science with ALFA. GALSPECT has a spectral resolu- 
tion of 0.18 km s _1 , and a fixed bandwidth of 1380 km 

5 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 
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s . Each of the 7 beams of ALFA has a 3.35 arcminute 
half-power beam width with a beam ellipticity of 0.2. 
The region was initially selected for significant structures 
in the local, low velocity 21-cm ISM, and the subsequent 
imaging of HVC 160.7-44.8-333 (hereafter HVC-A), HVC 
162.8-46.3-303 (HVC-B) and HVC 156.3-45.1-304 (HVC- 
C) was purely coincidental. The region was observed in 
a 'basketweave' or meridian-nodding mode, interlacing 
scans from day to day, sim ilar to observations detailed in 
IStanimirovic et alJ <|2006|) . 
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3. DATA REDUCTION AND ANALYSIS 

The raw spectra were calibrated for the IF bandpass 
using a frequency switching technique that employs a 
least-squares solution to distinguish the RF spectrum 
from the IF bandpass ( Heiles 2005a) . The spectra were 
then cleaned for baseline ripple caused by reflections in 
the optical fiber between the receivers and the spec- 
trometer. After Doppler-correcting the data to the LSR 
fra me, the data we re then corrected for fixed-pattern rip- 
ple (Hcilesi l2005bj) . Due to reflections in the Gregorian 
dome of the Arecibo telescope, ripples in the spectra of 
order 0.2 K amplitude can occur with a typical 'wave- 
length' of 300 km s~L This ripple is rather constant 
over the course of a few hours in this observing mode for 
a single beam, but varies strongly amongst beams. This 
ripple was mostly eliminated by fitting out residual ripple 
as compared to the overall average of each day's obser- 
vations. Each scan's data were then compared with all 
other scans' data at the sky position where they cross to 
determine the relative amplitude gains of each beam on 
each day. The data were re-calibrated for these relative 
gains, re-reduced for fixed -pattern ripple and then g rid- 
ded in a 2' x 2' map fsee IStanimirovic et alJ 1)20061) for 
details). The gridded data were then scaled to the equiv- 
alent region in the Leiden-Dwinge loo Survey (LPS) for 
a sin gle, overall gain calibration ijHartmann fc BurtonI 
1997). 

Once this final data cube was attained, the HVC was 
fitted with a single Gaussian at each pixel to determine 
basic characteristics of the cloud. In our data cube, mul- 
tiple velocity components do not occur in the same pixel, 
and the single Gaussian fit accurately reports the fluxes, 
velocities and line widths in the data cube. 

4. LHE HVCS 

The HVCs imaged are some of the highest Galactic 
standard of rest (GSR) velocity clouds known, at —280 
km s _1 , and can be referred to as v ery high velocity 
clouds (VHVCs). iPutman etaD <|2002|) showed that the 
distribution of HVCs in the HIPASS data set (all S < 2°, 
resolution of 15.5') can be modeled with a Gaussian dis- 
tribution centered at —38 km s -1 GSR with a dispersion 
of 115 km s _1 . The updated Wa kker fc van Woerdenl 
(1991) catalog (all sky, resolution < 30') has a mean 
of —45 km s^ 1 with a dispersion of 104 km s _1 .Thc 
catalog has been updated by i ncluding cloud s from the 
catalog of HVC components bv lMorras et alJ 1)20001) for 
5 < 23° and includes clouds with |VLSR| > 90 km s" 1 
Csee lWakkerl ll200l for more information on the catalog 
properties). It is clear that an HVC with a GSR ve- 
locity of —280 km s _1 is at an extreme of the velocity 
distribution, though not inconsistent with these samples. 
We observed 3 large, separated clouds, HVC-A, HVC-B 
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Fig. 1. — The image shows the HI column density and the cen- 
tral velocity along the line of sight, derived by modeling each HI 
spectrum with a single Gaussian function. Color represents the 
center velocity of the Gaussian in the GSR frame and brightness 
represents column density, as shown in the color bar (top). Pic- 
tured in full are three clouds, HVC 160.7-44.8-333 (HVC-A), HVC 
156.3-45.1-304 (HVC-C) and HVC 162.8-46.3-303 (HVC-B), and 
the incomplete edge of HVC 166.9-43.2-275 (HVC-X). Also pic- 
tured is the very faint HVC 160.0-48.4-286 (HVC-U). Features are 
also labeled that are used in our modeling. 



and HVC-C that have all been previously cataloged (see 
Figure ^) . These three clouds have a mean full- width 
at half-maximum (FWHM) for a single gaussian fit of 
28 km s" 1 , with a dispersion of 8 km s _1 . This is a 
typical FWHM value for HVC s, in agreement wit h the 
34 km s _1 median value from Ide Heii et all 120021 We 
do not have a complete image of the fourth large cloud 
that is visible in the LDS map of this sub-complex, HVC 
166.9-43.2-275 (HVC-X), so we will leave it out of the 
analysis. Also pictured is the previously unknown HVC 
160.0-48.4-286 (hereafter HVC-U). HVC-U is angularly 
compact (~ 30') and in any given pixel approaches our 
detection limit (~ 10 19 cm~ 2 ). Note that this is an ex- 
ample of the detection of a new HVC with ALFA, and 
confirms the claim that ALFA will detect new low col- 
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u mn density clouds in the Galactic halo 6 . 

iWeiner et alJ (|20D If) have inferred a distance to the 
observed sub-complex of 13+g kpc. This distance was 
measured using the Ha flux, under the assumption 
that all of this flux comes from ionization of HI due 
to the Galactic interstellar radiation field (ISRF) and 
subsequent radiative recombination. It is important 
to note, however, that the model for the UV flux in 
the halo does not take into account possible patchiness 
in this radiation field that could arise from variabil- 
ity in the underlying structure of sources of UV radi- 
ation (O and B stars) , self-shadowing of HVCs to reduce 
Ha flux or possible effects of colli sional ioniozation (see 
iBland-Hawthorn fc Putmanll2f)0ll for an in-depth discus- 
sion). This model also fails to predict the observed Ha 
brightness of the Magellanic Stream. In light of these 
possible sources of error, an independent confirmation of 
this distance result would significantly bolster the find- 
ing. 

Morphologically, HVC-A is distinctly elongated, with 
an aspect ratio of around 6. It also has a few small 'arms' 
(labeled f-1, f-2, and f-3 in Figure^) that are folded back 
along the direction of elongation of the cloud. These 
'arms' typically seem to show some amount of slowing 
as compared to the bulk of the cloud, as if they are 
being swept back by the motion of the cloud through 
the surrounding medium. This velocity gradient may 
be related to the action that generates 'head-tail' mor- 
phology observed in other HVCs, whic h may be caused 
by ra m pressure and viscous stripping (jQuilis fc Moord 
120011) . Each of these 'arms' points within the same 45° 
cone in the plane of the sky, strongly suggesting that 
the velocity component of HVC-A in the plane of the 
sky is toward the upper left of the figure, generally along 
the length of HVC-A. This is consistent with the fact 
that the surface of HVC-A is much smoother on this 
'u pwind' side, as if sculp ted like a raindrop (see models 
in lOuilis fc Moore) 1200 ID . We extrapolate further that 
the smaller three clouds were once part of HVC-A and 
that the action that is affecting the 'arms' on HVC-A 
has a similar effect upon HVC-B, HVC-C and HVC-U; 
all of these clouds have relatively similar velocities (see 
Table ^1 and are 'downwind' of HVC-A. We also detect 
a faint bridge from f-1 to HVC-B, which lends credence 
to these claims. Under these assumptions, we can learn 
much about the clouds' physical attributes through the 
study of gas drag. 

5. HVC DRAG 

To make headway in understanding the interaction of 
this cloud with its environment we need to make a num- 
ber of simplifying assumptions. It is not the authors' 
intent that these assumptions be exactly correct, but 
rather that they are plausible; in addition, they give us an 
idea of the accuracy of the result and the ways in which 
it is susceptible to systematic errors. Before we begin, let 
us define a 'feature' as any part of the cloud whose move- 
ment with respect to HVC-A we wish to model, whether 
it be a swept-back arm of the main body (i.e. f-1, f-2 
and f-3 in Figure 0) or a distinct cloud, such as HVC-B, 
HVC-C and HVC-U. Three key assumptions regarding 

6 see the GALFA whitepaper, 

http : / / www . naic . edu /alf a/ galf a/ docs/ galf a_white_sept25 . pdf 



the interaction of HVC features with the ambient halo 
gas to keep in mind in the following analysis are (1) that 
the surrounding halo gas that generates the drag on these 
HVCs is static enough that its motion is negligible, (2) 
that the variation in the direction of the 'apparent wind' 
that the features experience is small enough to be ig- 
nored, and (3) that the displacement of a feature from 
it s origin is governed by the effect of drag. 

Bcniamin & Danlv (1997) (hereafter BD97) proposed 
that the distances to HVCs could be determined by mak- 
ing the assumption that they were at terminal velocity. A 
cloud moving at terminal velocity in a gravitational field 
g can be described by a simple force-balance equation: 

M c g = ^C D p h v c 2 A, (1) 

where v c is the velocity of the cloud through the medium, 
Ph is the density of the medium it is moving through, and 
A is the area that the feature presents to the medium per- 
pendicular to the direction of motion (BD97). M c is the 
mass of the feature and g is the acceleration felt by the 
cloud due to Galactic gravity. Cd here is the 'drag coef- 
ficient', a way of parameterizing the dependence of drag 
forces on shape - aerodynamic shapes can have Cd < 0.1 
and bluff shapes can have Cd > 2. The basic premise 
of BD97 is that given (1) a fiducial model of the gas- 
density of the halo, (2) the assumption that all HVCs 
are at their terminal velocities, and (3) that their space- 
velocities are roughly equivalent to their observed veloc- 
ities, we can find distances to these HVCs. Applying 
this reasoning and the ISM models used in BD97 to the 
3 clouds of interest here, would put the sub-complex at 
a distance of ~4 kpc (see Table ^ and Figure [21 • Al- 
ternatively, if we assume the Ha-measured distance of 
13 kpc, the predicted density of the halo medium is an 
order of magnitude larger than the model assumed in 
BD97. It is possible to reconcile the terminal-velocity 
picture with the specific information pertaining to this 
HVC only by strong fine-tuning of the relevant parame- 
ters (e.g. setting the space velocity of the clouds to the 
Galactic escape speed and setting Cd = 2). 

In the following analysis, we do not assume that these 
HVCs are at their terminal velocities, though they expe- 
rience some drag from the Galactic halo. The equation 
of motion for a cloud in the Galactic halo that is subject 
to drag by the gaseous halo is 

M c a = ^C D p h v c 2 A + M c g, (2) 

where a is the acceleration. We have assumed that mag- 
netic fields are weak and do not significantly affect the 
motions of the clouds. Two clouds at approximately the 
same location suffer a differential acceleration 

A« = \c DPh A (£) * ICWA (£) , (3) 

where E c = M c /A is the surface density of the cloud, and 
where we have used the fact that the clouds have similar 
velocities. Clouds suffering a differential acceleration Aa 
will develop a differential velocity Av and a separation s 
related by 
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This assumes that the clouds started at the same velocity 
before separation. It follows from Eqns. |3and0]that the 
number density in the halo (measured in hydrogen atoms 



3 ) is 
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where iVn is the cloud column in hydrogen atoms cm~ 2 . 

To put this equation in terms of observables, we must 
introduce the angle formed by the overall cloud complex 
velocity vector with respect to the line of sight. We call 
this angle <j>. The observed velocity is then 



Kbsl = v c cosct>, 



(6) 



where w bs is negative in this case. The spatial separation 
of the clouds s is given in terms of their observed angular 
separation by 

0D = ssin<?!>, (7) 

where D is the distance to the clouds. Equation (J5J then 
becomes 
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In the definition of 2 the only two variables not mea- 
sured by observation are Co and <fi. We make the as- 
sumption of Co = 1 (as in BD97) as a fiducial estimate. 
To estimate <j> we take advantage of the fact that the HVC 
complex under investigation is at a very high velocity and 
so will typically have a smaller than the general HVC 
population; were the <p to be large the true velocity would 
be anomalously large (see eqn EJ and we would expect 
to see other HVCs with much larger observed velocities. 
To quantify this assertion we constructed a mock pop- 
ulation of HVCs that correspo nd to the velocity distri- 
bution observed in the updated Wakker fc van Woerdenl 
l)1991|) catalog. We give these HVCs a Maxwellian dis- 
tribution of true velocities so that when a naive random 
distribution of cj) is taken into account the mock velocity 
distribution reflects the broadly Gaussian catalog veloc- 
ity distribution. In the catalog, the HVCs in question are 
on the tail of the distribution, at a — 2.3. We investigate 
this same region in the mock population, and find that 
clouds at this same velocity have = 25 ± 10 degrees. 
We take 25 degrees as our fiducial value for </>, and note 
that a systematic error in <f> can result in errors in 5 of 
~ 40%. 

To reproduce the Gaussian distribution of observed ve- 
locities, this analysis makes the assumption that the true 
velocity unit vectors can point anywhere in space with 
the bias towards net infall. While this velocity distri- 
bution is wh at one would expect from models like those 
described in iMaller &: Bullock! l)2004|L these models are 
still far from general consensus. For rougher limits that 
do not depend upon models we can assume that the true 
velocity of this cloud does not exceed the highest veloc- 
ity ev er observed in an HVC in the halo, currently -450 
km/s (de Hcii ct all 1200^1 . This gives us a maximum <p 
of 51°. For </> less than this maximum value, there is only 
an 8% chance that <f> < 15°, our minimum value from the 
mock distribution. Even with this extreme range of 0, 2 
only varies by a factor of 2. 

To do our modeling, we apply this equation to six dis- 
tinct features that we believe may have been pulled off of 



HVC-A (see Figure [TJ. Three of these features, f-1, f-2 
and f-3, are small 'arms' connected to this main body, 
and certainly have the swept-back appearance that led 
us to this reasoning in the first place. The other three re- 
gions, HVC-B, HVC-C and HVC-U, are the three other 
distinct clouds in the map assumed to have been once 
connected to HVC-A. We neglect modeling HVC-X as it 
is not completely imaged. We do not attempt to model 
other features as difficulty in the evaluation of model pa- 
rameters (0, Aw b s or iVu) precludes their accurate anal- 
ysis. In particular, features with confusing morphology 
have no reliable value for 0. 

To determine a value for AA^ 1 we first note that our 
measurement is of the neutral hydrogen column density, 
rather than the total proton column density, so we use 
a fiducial correction fa ctor 4/3 to account for a small 
fraction of ionized gas l)Tuftdl2004|) . We then compare 
the average column density in each feature to the average 
column density in the 'head' of HVC-A, located at / = 
162.7, b — —43.2. The column density of the front of 
HVC-A, which is presumed to generate the drag, is the 
value of interest. Over the disruption timescale, 10 Myr, 
the cloud travels only ^5 kpc, so our implicit assumption 
that the density of the ambient medium is constant over 
the disruption time is not wholly unreasonable. 

The errors in the analysis are dominated by systemat- 
ics such as our measurement of the distance between the 
features and the main cloud and the unknown 3D orien- 
tation and movement of HVC-A and the other features. 
Since the values of 2 we predict vary with a fractional 
standard deviation of ~ 0.4, we take this as a reasonable 
value for our random errors, of the same order as our 
predicted errors in <p and significantly smaller than the 
most conservative errors in <j>. We find the median value 
to be 2 mo( j = 2.7 cm -3 pc. The results of this analysis 
are listed in Tabled 

6. DENSITY ALONG THE LINE OF SIGHT 

To determine the validity of these results we must 
compare them not only to existing distance estimates, 
but also to estimates for the density of the ISM along 
the line of sight to the HVCs. In the Galactic disk 
regime we simply assume a plane-parallel Cold Neutral 
M edium (CNM) and Warm Ioni zed Medium (W IM), as 
in lDickev fc Lockmanl (|199(T l and lRevnoldsl l|1993lL Since 
this region is observed at b — —45°, once the plane- 
parallel assumption becomes questionable along the line 
of sight we will have left the disk regime and must con- 
cern ourselves only with the Hot Ionized Medium of the 
halo (HIM), which dominates the density. 

In BD97 the assumption is made that the model for 
the halo at R = Rq determined in W95 can be sim- 
ply extrapolated as a plane-parallel medium, as we have 
for the other phases of the ISM. In our case, we can- 
not make this simple extrapolation for the halo, as our 
object lies at R ~ 20 kpc, where the character of the 
halo is rath er different than it is locally. In a sepa- 
rate paper l)Peek et alJ 12006ft . the authors incorporate 
result s from OVII absorption measur ements llFang et al. 
2006), halo emissivity measurements l|Kuntz fc Snowderi 
2000), and recent th eoretical constraints on disk pressure 
(Wolfi re"et al]l2003Tl to develop a self-consistent, isother- 
mal hydrostatic halo model. The model is essentially the 
simplest possible model of an isothermal gas in a known 



HVC Distances Through Morphology 



5 



10° 





istance 








H alpha d 














|_|| | Djc ^ ej , & Lockman 1990) \ \ 

Warm HII (Reynolds 1993) \ 

Concordant Halo (Peek el al 2006) • 






\ \ * 


Total ISM density i 








10 2 10 3 


10 4 


10 5 



D(pc) 

Fig. 2. — This Figure is a plot of total ISM hydrogen num- 
ber density, njj vs. distance towards the HVC complex. The 
theoretical density profile is a sum of the neutral HI disk from 
Dicke yfc Lockmanl ^990), the warm HII Reynolds layer from 
Reynolds! ( 119931) and the concordant halo model from Peck ct al. 
2006, in prep. The horizontal lines in the lightest shaded area are 
the predictions for the density of the ambient medium from the 
assumption of terminal velocity. The vertical line is the distance 
as predicted by Ha mea surements, with the associated shaded area 
representing error bars (Wciner ct al. 2001). Our dynamical esti- 
mate, S me< j is plotted in a dashed line; the dark shaded region 
represents estimated errors of a factor of 2. The asterisk in the top 
left of the plot represents H max (see Equation |5J. Note that the 
intersection of the modeled nn(D) with the measured distance is 
consistent with our result, and that the terminal velocity distance 
is not. 



gravitational potential that is consistent with modern ob- 
servational results. It is this global halo model to which 
we compare the Ha distance and our dynamical results 
in Figure |21 

7. DISCUSSION 

The results for S are in Tabled along with the average 
and peak column density of each feature. These values 
vary by a factor of few, which is not surprising consid- 
ering the assumptions with regard to morphological de- 
tails of these clouds along the line of sight. In Figure 
[21 we show the halo number density of hydrogen nn as 
a function of dist ance towards t he clo ud, as well as the 
Ha distance from lWeiner et all 1)200 lj) . We also plot the 
prediction, S, from the dynamical analysis, along with 
the standard deviation within our sample. Our predic- 
tions are in good agreement with the intersection of the 
cumulative density profile and the Ha distance. 

In future observations, we can use our measured quan- 
tity, S, to give us a handle on the distance to a cloud, 
if we assume the canonical density models in Figure |2 
Below about 2 kpc, where we take the halo to be domi- 
nated by the warm Reynolds Layer and the HI disk, the 
models show n# oc Z3~ 3 / 2 , so that S = nnD oc D~ x / 2 
and our prediction for distance would be D oc S -2 . This 
implies that fractional errors in our quantity will 
result in fractional errors in the distance 5D/D oc 25H/3. 
Evidently, our capacity to estimate distance is rather 
hampered in this nearby regime, as the errors that arise 
from confusion and projection effects are strongly ampli- 



fied. This issue is of little relevance, though, as there 
are few direc t measurements of HVCs so nearby (see 
lWakkerll200ll) . The picture in the distant regime, 2 
kpc < D < 100 kpc, is much rosier. Since in this 
region, hh oc D~ 5 / 2 , we have D oc S~ 2 / 3 , so that 
SD/D oc 2(5S/3S. This means that our errors in this fig- 
ure lead to reduced errors in our distance measure, which 
gives a powerful measure of distance. Alternatively, for 
the few clouds that have good distance measurements, 
this technique could be turned around to measure densi- 
ties and test theoretical models of the halo density profile 
and smoothness. 

A concern with this kind of technique is that it could 
be applied to HVCs coincidentally along the same line 
of sight and thereby generate incorrect distance (or den- 
sity) estimates. This is of course possible, but if the 
clouds have dissimilar velocities, S will increase dramat- 
ically. As 5 rises, the predicted distance decreases, but 
as the Galactic disk is approached the density profile, 
dominated by the HI disk, flattens out and the dynam- 
ical prediction no longer intersects it in any place. This 
defines the maximum possible value of 2: 

„ 39 cm -3 pc 

"max — 7 > \") 

COS 

for the assumed IDickev fc Lockmanl (1990) HI disk den- 
sity profile (see asterisk in Figure [2J). Therefore, we can 
easily get an idea if it is impossible for two clouds to be 
related by drag forces - if their measured S is greater 
than this maximum value, then they cannot have been 
sheared apart by drag. 

We have shown that modeling the movement of various 
pieces of our HVC complex as having been stripped from 
the main body of the cloud by simple drag forces is con- 
sistent with reasonable halo models and the inferred Ha 
distance. One exciting aspect of this conclusion is that 
it bolsters the Ha distance measure technique. Given 
both the large uncertainties regarding the magnitude and 
patchiness of the radiation field, and some uncertainty in 
whether the Ha is generated solely by photoionization or 
whether collisional ionization plays a role, it is important 
that there exist checks on this distance measure. Since 
field patchiness can be stronger closer to the disk, our 
confirmation of the Ha distance technique serves to show 
that patchiness is a limited effect at z ~10 kpc or higher , 
in agreement with iBland-Hawthorn fc Putmanl l)2001[) . 
at least toward this cloud. Collisional ionization is ex- 
pected to increase with HVC velocity and halo density. 
Given the extreme velocities of these clouds, we claim 
that collisional ionization does not strongly effect the Ha 
distance measurement for any HVCs in a similar density 
regime. 

We have inferred a distance above the plane for these 
clouds of ~ lOkpc. This calls into doubt whether these 
clouds could have been produced by a Galactic fountain 
model, where gas blow-out from the disk HIM launches 
material in to the halo, whi c h later reconden ses and falls 
as HVCs (jBregmanl Il980l Ide Avilled 1200(1) . Galactic 
fountain models set a maximum height for upflowing gas 
that condenses into HVCs at <10 kpc. Therefore any 
clouds at this height should be moving relatively slowly, 
as they are at the apex of their cycle, which contradicts 
the large velocities we measure. Also, at this distance, 



() 



the complex has a total mass of 2.7 x 10 5 M© and a kinetic 
energy of 7.2 x 10 53 ergs. This is equivalent to the en- 
ergy output of approximately 720 supernova explosions 
into the ISM. 

At our inferred position, W95 predicts a gravitational 
potential of —1.16 x 10 5 (km s -1 ) 2 . This implies that if 
the cloud started at rest far from the Galaxy and con- 
served energy as it fell, it would have a velocity of 480 
km s _1 . As this velocity is in excess of the estimated 
space-velocity of the HVCs, we find that our distance is 
consistent with models in which HVC velocity stems from 
cloud free-fall in the Galactic potential. If this particular 
group of HVCs started from rest we can also determine 
that the HVCs could not have originated any closer than 
~ 60 kpc. 

The inferred density at this position is certainly in 
agreement with the observational constraints discussion 
in Section H3 as they are used to construct the halo 
model to which we compare our results. The density 
is also in rough agreement with inferred halo densities , 
such as those put forwar d bvlWeiner fe Williams! |l996), 
n h > 10- 4 cm-3, and iQuilis fc Moord l|2001|) . n h (10 
kpc) < 3 x 10~ 3 cm~ 3 . 

Models in which HVCs form f r om the cooling of halo 
gas (e.g. IMaller fc Bullockll2004l IKaufmann et alJl2006L 
Sommcr-Larscn 2006) allow comparison both in the HVC 
parameters and in the halo density. The HV C dis- 
tance scale determined in IKaufmann et all ||20O6) of 10- 
20 kpc is in line with o ur observations and simulations in 
iSommer-Larsenl l)2006l) have cloud masses consistent with 
our derived cloud masses. Sommcr -Larsenl ll2006|) also 
finds hot halo densities 10 -3 5 cm -3 > uh > 10 -4 cm" 3 
at R = 20 kpc, and IMaller fc Bullockl (|2004ft find n H = 
2 x 10~ 4 at R = 20 kpc in their "hot gas core" , both in 
good agreement with our results. It seems that our mea- 
surements are in good agreement with this theoretical 
picture of HVCs and the halo. 

8. CONCLUSION 

We have demonstrated that we see the effects of halo 
drag on a group of HVCs in their morphological and kine- 
matic structure. Indeed, a simple model of drag physics 
gives a consistent and plausible result for the distance to 
these HVCs. This distance is in agreement with the Ha 
distances measured for these clouds. At this distance, 
we confirm a halo density of njj — 2 x 10~ 4 cm~ 3 at 
R ~ 20 kpc. These results are in rather good agreement 
with the modern theoretical picture of a circum-Galactic 
origin for HVCs, and in rather strong disagreement with 
typical Galactic fountain models. 

To reproduce this analysis in other HVCs, one would 
need to both be able to identify similar morphological 
clues in the structure of HVCs (to determine 0) and be 
able to estimate </> with similar precision. While this is 
certainly possible in some other HVCs, we predict that 
the vast majority of HVCs will elude this technique in its 
current form. Whether the physics described herein in- 
fluences the structure of larger HVC complexes and can 
be used to measure HVC distance and halo density has 
yet to be shown. In the few cases where this technique is 
usable, comparison to simulation would allow us to eval- 
uate whether more sophisticated analyses are required to 
model fully the physics involved, or whether HVC distor- 
tion and destruction, for typical physical parameters, can 



be robustly modeled with our naive gas-drag analysis. It 
is also interesting to note that the structure of a hydro- 
static halo is a strong function of temperature on 10 kpc 
scales so any future measurements of S may help tightly 
constrain halo temperatures, and subsequently models 
of spiral galaxy formation and evolution. We look for- 
ward to more in-depth quantitative analysis of cloud and 
cloud-complex morphology in the presence of the Galac- 
tic halo. 

The authors would like to acknowledge those who 
helped get GALFA going, including Jeff Mock, Aaron 
Parsons and Dan Werthimer as well as the Arecibo staff. 
JEGP would like to acknowledge many helpful conversa- 
tions with Evan Levine and Leo Blitz. The research of 
CFM is supported in part by NSF grant AST00-98365 
and NSF grant AST-060683. The research of CH and 
JEGP is supported in part by NSF grant AST04-06987. 



HVC Distances Through Morphology 



7 



TABLE 1 

Parameters for HVCs and features imaged. The parameters were determined by setting 
boundaries for the clouds where the error in the value for the height of the gaussian fit to the 
spectra reached 25%. the positions noted are the geometric centers of each of the features, and 
do not necessarily reflect the inferred values of 0. features f-l, f-2 and f-3 were also bounded 
by a cut in velocity to distinguish them from the hvc-a. this cut was chosen at a minimum of the 

velocity histogram in the region. 



Name 1 b VGSR (km s" 1 ) (Nh) (l0 19 cm~ 2 ) N H , pca k (l0 19 c m - 2 ) G(°) H (cm~ 3 pc) 



HVC-A 


160.7 


-44.8 


-280 


6.0 


15 






HVC-B 


162.8 


-46.3 


-257 


4.8 


9.1 


3 


2.8 


HVC-C 


156.3 


-45.1 


-240 


2.8 


4.1 


4 


2.3 


HVC-U 


160.0 


-48.4 


-236 


1.8 


3.0 


2 


3.0 


f-l 


161.6 


-44.9 


-270 


2.2 


4.8 


0.5 


0.80 


f-2 


160.7 


-43.4 


-268 


2.5 


7.2 


0.4 


1.8 


f-3 


157.3 


-46.6 


-267 


1.8 


3.7 


0.4 


1.3 
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